State-dependent impedance of a strongly coupled oscillator-qubit system 
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We investigate the measurements of two-state quantum systems (qubits) at finite temperatures 
using a resonant harmonic oscillator as a quantum probe. The reduced density matrix and oscillator 
correlators are calculated by a scheme combining numerical methods with an analytical perturbation 
theory. Correlators provide us information about the system impedance, which depends on the qubit 
state. We show in detail how this property can be exploited in the qubit measurement. 
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I. INTRODUCTION 

The growing interest in quantum information theory 
and the rapid development in nanotechnology have re- 
sulted in an extensive study of quantum two-state sys- 
tems (qubits) and the quantum measurement theory. 
The measurement theory of frequency independent de- 
tectors has been studied in detail in recent yearsiSii^ 
Motivated by the new directions of research, we investi- 
gate possibilities to use a harmonic oscillator as a generic 
frequency dependent measuring device of a qubit. 

Our system under study consists of a qubit coupled res- 
onantly to a harmonic oscillator (Fig. [I}, both coupled to 
a bosonic heat bath at a finite temperature. This model 
has a wide range of applications in solid-state physics as 
well as in quantum optics and has raised considerable at- 
tention ! a telyi£i§iL2i2ii2iiiii2 In the literature there exists 
various propositions to realize this system. In solid-state 
physics the harmonic oscillator is realized by a resonator 
circuit and the qubit, for example, by a Josephson charge 
or flux qubit. The heat bath corresponds to the elec- 
tromagnetic environment of the circuit. The connection 
of these system s to c avity QED has been explained and 
studied in Refs. loll.il 

Recently, also experimental studies of solid-state real- 
izations of cavity QED have become accessible. Wall- 
raff et a/^i performed an experiment where they suc- 
cessfully managed to couple a transmission line resonator 
with a Cooper pair box acting as a qubit. They found 
a clear evidence of a quantum entanglement between 
the resonator and the qubit by measuring the vacuum 
Rabi mode splitting. Chiorescu and coworkers^ realized 
the oscillator+qubit system by coupling a Josephson flux 
qubit to a SQUID. In their experiment the SQUID be- 
haved like a harmonic oscillator and was used as a mea- 
suring device of the qubit. The entangled states of the 
oscillator-qubit system could be generated and controlled 
through the use of microwave spectroscopy. The fact that 
the qubit is fixed on the same chip as the resonator al- 
lows circuit QED systems to explore the strong coupling 
regime of cavity QED, difficult to reach in quantum op- 
tics realizations^. 



In the light of these recent experiments, we examine 
how the qubit information can be extracted from the os- 
cillator measurement. Oscillator correlators reveal the 
entanglement between the qubit and the oscillator. The 
correlator information can be linked to experimentally 
accessible quantities. We propose a scheme to extract 
information about the state of the qubit by probing the 
impedance of the oscillator. This can be realized by a 
transmission measurement. As a result of such an exper- 
iment, one can resolve the diagonal state of the qubit. 
Such transmission measurements have been recently ex- 
ploited in studying the solid-state cavity QEDii. 

Our paper is organized as follows. In Sec. II we de- 
tail the general formalism for calculating the correlators 
and density matrices. We also explain the basics of the 
susceptibility measurement, and discuss the experimen- 
tal realization of this measurement. In Sec. Ill we study 
the susceptibility in various circumstances and how the 
state of the qubit can be read from the susceptibility 
information. Finally, we summarize and discuss the sus- 
ceptibility measurement in Sec. IV. 
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FIG. 1: Possible realization of the studied system. The res- 
onator circuit is coupled to the qubit represented by a. In 
practice a could be for example a Cooper-pair box with a 
capacitive coupling to the oscillator. 



II. GENERAL FORMALISM 

Our starting point is the Jaynes-Cummings Hamilto- 
nian, which describes a qubit coupled to a harmonic os- 
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cillator in the absence of dissipation: 
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+ ^a x (a + J). (1) 



Operator a/a operates on the qubit /oscillator degrees of 
freedom. The operators a and a' obey the usual bosonic 
commutation relations. We assume the oscillator in res- 
onance with the qubit, ujq = aj q b- In such a case the ex- 
citation spectrum is doubly degenerate if oscillator-qubit 
coupling g vanishes. A finite coupling g <C u>a slightly 
lifts the degeneracy of the spetrum leading to the Rabi 
splitting. The properties of the spectrum of Q are ana- 
lyzed in detail in Ref. IT^l 

The total Hamiltonian under study is 
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where = YnhuJibJbi and = Hjfiuij c^Ci de- 

scribe the environments of the qubit and the oscillator, 
respectively. In addition, iff, = t7 x Y*igi(bi + bj) and 
i? ; ° s t c = {a) + a)T,jgj(cj + c^) couple the system to the 
environment. The chosen bath model is a popular choice 
for its conceptual and technical simplicity^. 

Our strategy is to calculate the reduced density matrix 
of the oscillator-qubit system and the oscillator corre- 
lators by numerically diagonalizing the qubit+oscillator 
system and taking the environment into account pertur- 
batively. Only the lowest-order contributions of the per- 
turbation series are assumed to be significant. However, 
to obtain a finite linewidth of energy levels of the os- 
cillator+qubit system, the perturbation series must be 
analytically resummed to the infinite order. 



A. Perturbation theory 

We apply the perturbation theory for an open system, 
developed, for example, in Refs. IllmqllSl . The time evo- 
lution of the density matrix involves forward and back- 
ward propagators which get coupled after tracing out the 
environmental degrees of freedom. Given a complete set 
of energy eigenstates of a quantum system, the evolution 
of the density matrix can be expanded in this basis as 



(3) 



where TL n ^ n }(t,t') is the reduced propagator and 

fmn 2 {t') is the initial reduced density matrix. In the 
graphical language the propagator is the sum of all Feyn- 
man diagrams with the given states at each end. The 
Feynman rules for an open system, particularly for a 
quadratic environ ment w ith a bilinear coupling term are 
described in Refs. [Till 

Supposing that at the initial time the density matrix 
of the system+environment is of a tensor product form, 



the perturbation theory is significantly simplified. In the 
frequency space the propagator is defined as 
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IQ& t')d(t-t'). (4) 



The full propagator is described by the Dyson equation 

n nXn Uu;) = u onln Hu)+ 

712 Tin 712 n 2 

n^Ms^Mn^M, (5) 

where T,^(ui) is the irreducible self-energy and II (ui) 
is the free propagator. Summation over the repeated in- 
dices is implied. The free propagator elements are deter- 
mined by the energy eigenvalues of the reduced system 
alone: 



U° nln '}(tu) = 



(E m — E ni )/h - 



i 1 ^ 



In our case, the energies in Eq. 10 are the exact energy 
eigenvalues of the oscillator+qubit system. The pertur- 
bation theory is needed only in taking the external baths 
into account. After calculating to the desired or- 

der, the exact propagator can be solved from Eq. JSJ. 
In our work we only take into account the second-order 
contributions to S*?(w), see Fig. [21(b). For example, the 
first of the self-energy diagrams represents the expression 
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, u + (E n2 - E n )/h + ui' + ir) 
x (n\x\nl)(n3\x\n)6 n2n 4, 



where J(u/) is the spectral density of the bath and /(u/) 
is the Bose distribution function. The spectral density is 
chosen to produce an ohmic damping. 

Correlators (A(t 2 )B(ti)) contain important informa- 
tion about the system. Generally, correlators are func- 
tions of the two time variables t± and t 2 , but in a steady 
state they depend only on the relative time t 2 — t\ due 
to the temporal invariance. In this paper we make use of 
a mixed representation where correlators depend on one 
time variable t — min (t\, t 2 ) and one frequency variable 
lu: 



D AB (co,t) 



^(AfaMt^dlh-h]. (7) 



If the system is in a steady state, DAB(w,i) is indepen- 
dent of t. The reason for considering correlators of the 
type DAB(u,t) is that in Sec. IHII we study the temporal 
evolution of the quasistatic susceptibility of the oscilla- 
tor+qubit system. 

Correlators cannot be calculated from the reduced 
propagator and the initial density matrix alone. The di- 
agrammatic method of calculating correlators in an open 
system is presented in Refs. Il7lll8l To calculate the cor- 
relator (7J, we first need to solve the temporal evolution 
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of the reduced density matrix from the initial time t' to 
t = t\ (assuming t\ < t 2 ). This task can be performed 

with the help of II™ 1 "3/ (^i')- Next we have to calculate 
the vertex corrections in the first external vertex B, and 
propagate the system with II™ 1 "J,, (lj) to the other vertex 
A. Correlator (0 is schematically illustrated in Fig. 

In a stationary state one further simplification occurs. 
Quantum regression theorem^ allows one to calculate 
the correlator Q with the propagator and the initial den- 
sity matrix alone: 
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In the diagrammatic language this means that the ver- 
tex corrections in the external vertex B cancel out. When 
calculating the time-dependent correlator Q the vertex 
corrections must be included, see Fig. Since we cal- 
culated the self-energy to the second order, we take also 
the vertex corrections into account to the second order. 
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FIG. 2: (a) Propagator cvolVes-'the initial state specified by 
the density matrix f n p„2' ■ (b) Second-order contributions to 
the irreducible self-energy E. 



B. Numerical scheme 

For computing the propagator, the energy eigenvalues 
and eigenstates of the oscillator-qubit system are calcu- 
lated numerically. This is done by taking into account 
only the N lowest energy eigenstates of the oscillator. 
The lower the temperature with respect to Hloq, the less 
states N is needed. The propagator contains 167V 4 terms. 
Fortunately, the interaction with the bath couples only 
states of the reduced system differing roughly by Hloq in 
energy 1 ^. For this reason, the effective number of propa- 
gator elements is proportional to TV 2 . 

The reduced density matrix and correlators can now 
be computed from Eqs. and © (see also Fig. [2J ■ 




FIG. 3: Diagrammatic representation of the correlator £J. 
First the initial state develops to the moment t, then follows 
external vertex B and the vertex corrections. The last step is 
the propagation in the frequency space to the other external 
vertex A. 



C. Oscillator measurement 

In studying the effects of the qubit on the oscillator we 
need to identify convenient measurable quantities. The 
idea is to analyze how, by measuring oscillator observ- 
ables, we can extract information about the qubit. In a 
steady-state situation the cavity susceptibility of cavity 
QED is defined as 



'{[x{t),x{Q)])dt, 



(9) 



where x = (a + a'). The direct analogue of the x- 
susceptibility © in solid-state applications is the in- 
susceptibility of a resonator circuit. The resonator cir- 
cuit is characterized by two parameters, the inductance 
L and the capacitance C. Defining = y 'L/C and fol- 
lowing the conventions of Ref. |20| the (/(-susceptibility is 



where 



hZh 



[0(t)>(O)])di, 



(a + a r ). 



The conjugate of <f> is the charge operator 
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According to the linear response theory, the susceptibility 
is proportional to the impedance of the system, Z(lo) = 

To generalize to nonsteady-state cases, we also 

define a time-dependent susceptibility as 



X(w,t) 



e -(*'-*)([^'),^)])d(t'-t). (13) 



When relaxation caused by the environment is weak, 
x(ijJ,t) changes slowly compared with the internal dy- 
namics of the system and may be considered as a qua- 
sistatic quantity. This is the case with a high Q factor 
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system. In a nonsteady state the emission and absorp- 
tion between the system energy levels does not obey the 
detailed balance condition and x( w , t) may differ qualita- 
tively from the steady-state susceptibility studied in 
Ref.[l3|. However, the oscillator-qubit system approaches 
the thermal equilibrium irrespective of the initial state 
and, therefore, 



tion theory yields an approximate transition time 



lim x(^,t) = xH- 



(14) 



In practice, x(tL>,t) ~ x( w ) f° r t ~ Q/u>o- 

The susceptibility function contains knowledge about 
the dissipative and reactive response of the system un- 
der a weak external perturbation. The susceptibility at 
a certain moment x(u;,to) can be measured by expos- 
ing the system to a perturbation acting rapidly com- 
pared to the natural relaxation time of the system. In 
the solid-state context this could be carried out by con- 
necting the system of interest between input and output 
transmission lines (Fig. and measuring the transmis- 
sion of a microwave pulse. Different quantum states of 
the coupled oscillator and qubit system result in different 
impedances. By measuring the voltage of the transmit- 
ted signal one can acquire knowledge about the qubit 
state. The most important requirement for the single- 
shot measurement of this kind is that the mixing due 
to the microwave drive and the relaxation caused by the 
environment should take much longer than the voltage 
measurement necessary to identify the state of the sys- 
tem. Assuming the drive-induced mixing to be signifi- 
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FIG. 4: Realization to carry out transmission impedance 
measurement of a solid-state quantum circuit. The impedance 
Z(oj) represents the coupled oscillator-qubit system. The 
impedance Z(uj) is dependent on the quantum state of the 
system which can be probed by a microwave pulse. 

cantly quicker than the environment-induced relaxation, 
one can calculate a simple estimate for the mixing time. 
The driving Hamiltonian is assumed to take the form 



H d = V qcos{u d t), 



(15) 



where q is the charge operator ((12(1 and Vq describes the 
strength of the voltage induced by the microwave held. 
We assume that the frequency of the measurement signal 
tod coincides closely with an energy difference of two ar- 
bitrary energy eigenstates \i) and |/) of the system. An 
application of the first order time-dependent perturba- 
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(16) 



where Rq — h/e 2 w 25.6kf2 andp = i(a'—a). After time 
Td the microwave pulse has disturbed the system and lost 
the information about the initial state in xi^d, to)- 

The measurement time T m is limited from below by 
the voltage noise of the output signal. The measurement 
time can be estimated by the signal-to-noise formula 

AV 

V by-lm 

where Sy is the spectral density of the voltage fluctua- 
tions and AV is the voltage difference in the signal caused 
by two differing impedances Z\{lo) and Z 2 {u}) of the os- 
cillator. The different impedances Z\(iS) and Z 2 (u>) cor- 
respond to two different quantum states of the oscillator- 
qubit system. Supposing that the dominant voltage 
fluctuations originate from the amplifier, we can write 
Sy = 2k B T^ Zq, where Tjv is the noise temperature of 
the amplifier and Zq is the characteristic impedance of 
the transmission line. Setting S/N = 1 and solving T m 
from Eq. fT7| we get 



2fc B Tjv Z hk B T N Z 

T m = , . = 47T- 



(Avy 



{eAVYR Q - 



(18) 



The transmitted voltage difference between two different 
states is 



AV=\V 1 -V 2 \=V z, 



where we defined 



2(Z 1 /Z - Z 2 /Z ) 



(2Z 1 /Z Q + 1)(2Z 2 /Z + 1) 



(19) 



(20) 



The voltage difference is significant when Z\ and Zi are 
of order Zq or one is much larger and the other much 
smaller than Zq. 

For carrying out a transmission measurement capable 
of resolving two quantum states the condition Td T m 
must be fulfilled. Using Eqs. 1)16(1 and 1)18(1 their ratio is 



hV z 2 
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(21) 



In the analogous case of current driving 21 , the ratio be- 
tween the relaxation and measurement times is 



T m 



hI z 2 Z 



Zq 



2^z 2 



(22) 



k B T N \{Mi)\ k B T N JZ^ \(f\x\i), 

where x = + a and Iq is the amplitude of the driv- 
ing field. The matrix element in Eq. ((21(1 and 
in Eq. 1(22(1 are of the order of unity in transi- 
tions between the lowest energy states. The impedance 
Zq can be matched with Z\ or Z 2 via a matching circuit, 
in order to produce an optimal signal. If the requirement 
Td T m is not met, one needs to perform a sequence of 
measurements to resolve the qubit state. 
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FIG. 5: Real part of impedance at T = huio/k^ (red), T = 
ftwo/3fcb (blue) and T = ftajo/lOfcb (black). 



III. RESULTS 

Let us now turn to the impedance Z(u>, t) = i toxi^i t) 
of the strongly coupled qubit-oscillator system. We char- 
acterize the different qubit states by using impedance 
information and apply this knowledge to the qubit mea- 
surements. 

We illustrate the dependence of the oscillator 
impedance on the qubit state by considering an exam- 
ple. In the following numerical results, we choose the 
coupling strength g = 0.03wo between the oscillator and 
the qubit and the quality factor Q = 10 4 of the oscillator. 
These parameter values are accessible in the solid-state 
realizations^. The ratio between the mixing times 7 and 
k in the qubit and the oscillator, respectively, is chosen 
to be j/k — 0.08. This in agreement with Ref. 

The real and imaginary parts of the impedance Z of 
the system are shown in Figs. |5] and for three different 
temperatures T of the bath. The coupling between the 
oscillator and the qubit results in multiple peaks in strong 
contrast to a single oscillator. Moreover, the impedance 
of the coupled system is now strongly dependent on tem- 
perature. In the low temperature limit there are only two 
peaks in Re[Z(u;)], corresponding to the vacuum Rabi 
splitting. They have been experimentally observed re- 
cently for example in Ref. [n[ The susceptibility is stud- 
ied as a function of temperature in detail in Ref. IT^ . 

When the oscillator-qubit system is not in thermal 
equilibrium, the susceptibility changes radically (see 
Fig. 0). The detailed-balance condition, which relates 
emission and absorption processes, is violated. For ex- 
ample, one can have net emission of energy in some fre- 
quencies in addition to the absorption. In those frequen- 
cies the susceptibility takes negative values which lead to 
negative impedance peaks as in Fig. \ 7\ h) ■ The negative 
peaks are signs of spontaneous emission of energy and 
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FIG. 6: Imaginary part of the equilibrium impedance at T — 
huo/kb (red), T = hu>o/3kb (blue) and T = ?kJo/10fcb (black). 



relaxation towards the equilibrium state. 

Both sets of impedance curves in Fig. (a) and (b) 
contain characteristic peaks which can be utilized to re- 
solve the state of the qubit. At low temperatures the 
peaks III and IV are absent in Fig. (a). In addition, 
peaks I and II are negative in Fig. (b) . The transitions 
I-IV are shown in Fig. [S] They correspond to the tran- 
sitions between the lowest energy states of the coupled 
system. When the coupling g between the oscillator and 
the qubit is weak, the lowest energy states are given in 
terms of the uncoupled eigenstates as 
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The reason for the negative peaks in Fig. [7| (b) is that 
the ground state is unpopulated in the initial state and 
the system begins to spontaneously increase the ground 
state population. All peaks that remain present even in 
the zero temperature limit are higher and narrower at 
lower temperatures. At higher temperatures new peaks 
arise. Comparing curves Fig. 0(a) and Fig. 0(b) one can 
also notice that the growth of temperature decreases the 
differences between the impedance curves corresponding 
to the different states of the qubit. This phenomenon 
restricts the accuracy of the qubit measurement based 
on the impedance difference and puts an upper limit to 
the measurement temperature to roughly T ~ hug/kb. 

In Figs. and we have plotted the difference in 
the transmitted voltage amplitudes (|19|l . correspond- 
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FIG. 8: Transitions corresponding to the four peaks of the 
low temperature impedance curves (blue) in Fig. |7[a) and 
(b). The transitions I and II correspond to the vacuum Rabi 
splitting in (a). The reason that the peaks I and II are nega- 
tive in Fig. E^b) is that the initial state is far from equilibrium 
and spontaneously emits energy. 
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FIG. 7: (a) Real part of Z(u), to) at temperatures T — 
Tuvo/Wk b (blue), T = hu>o/2k b (red) and T = Tiuio/kh (green). 
The initial state at t = to is prepared so that the oscillator 
is in the thermal state and the qubit state is up (the lower 
energy qubit state), (b) Same as (a) but the initial state is 
prepared so that the qubit state is down. 
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FIG. 9: Difference between the transmitted current ampli- 
tudes corresponding to the qubit states up and down when the 
oscillator is in thermal equilibrium at T = huo/lOkh (blue), 
T = huo/Akb (red) and T = Tiuo/2k b (green). The ratio of 
the impedances Zt, and Zo is set to Zo/Z^ = 2. 



ing to the current (^-susceptibility) and voltage (q- 
susceptibility) driving, respectively. In Fig. the maxi- 
mum difference do not coincide with the peaks I-IV. This 
is because the impedance of the system at the peaks 
is much larger than the characteristic impedance of the 
transmission line Zq and the system effectively decouples 
from it. However, using an appropriate matching cir- 
cuit the effective impedance of the system can be varied. 
Better matching is obtained in Fig. llOl where we consider 
voltage driven system and the impedance ratio is chosen 
to Zq/Zi, = 1/500. The transmitted voltage amplitude 
can exceed unity because of the emission associated with 
the transition I and II. At temperatures T — huJo/lOk\, 
and T = tkoo / 4fcb the absolute values of the impedances 
at the peaks I and II are so high that the system again de- 



couples from the transmission line. Thats why the trans- 
mitted voltage corresponding to T = hLOo/2k^ gives the 
maximum for this particular impedance matching. The 
curves are not quite symmetrical with respect to loq be- 
cause the imaginary part of the impedance is antisym- 
metric with respect to ojq. 

The impedance of a slightly off-resonant oscillator- 
qubit system is plotted in Fig. 1111 The vacuum Rabi 
splitting, corresponding to the transitions I and II, is 
sensitive even to a slight detuning and thus the positions 
of these peaks follow the frequency oj q b of the qubit. On 
the other hand, the positions of the peaks corresponding 
to the transitions III and IV follow rather the frequency 
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FIG. 10: Same as Fig. but with voltage driving 
(impedances replaced by admittances, and calculated from 
the (/-susceptibility) and Zo/Zb — 1/500. 
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FIG. 11: Equilibrium impedance of a resonant and a slightly 
off-resonant oscillator-qubit system at T = 7hiio/2kb- The 
blue curve represents the exactly resonant case, the red curve 
corresponds to the case w q b = 1.01 too and the green curve to 
the case tj q b = 0.99 cjq- 



u>o of the oscillator. With small detuning, the height of 
the impedance peaks is only slightly altered. 

As the system relaxes towards the equilibrium state, 
the impedance settles to the equilibrium pattern. In 
Fig. E| we plot the temporal evolution of two nonequilib- 
rium impedances. They correspond to the initial states 
where the oscillator is in thermal cqulibrium and the 
qubit is prepared either up or down. The susceptibility 
changes slowly compared to loq 1 and can be considered as 
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FIG. 12: Relaxation of the impedance at T — ?ia>o/L8&b as 
a function of time. The initial state at t — is prepared so 
that the oscillator is in the thermal state and the qubit state 
is up (red line) or down (blue line). The figures correspond 
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quasistatic. After the time t > Q/uJo the susceptibilities 
in both cases are nearly equal, reflecting the uncertainty 
about the state of the qubit (Fig. El- When 7<k, the 
oscillator dissipation yields the dominant time scale for 
the relaxation of the qubit near resonance. 

IV. CONCLUSIONS 

In this paper we discuss the susceptibility and the 
impedance of a resonant oscillator-qubit system at finite 
temperatures. We have studied the response of the sys- 
tem in various noncquilibrium states and explored possi- 
bilities to apply this phenomenon to determine the state 
of the qubit. In certain conditions one can carry out a 
transmission measurement capable of resolving the qubit 
state. 

Let us estimate the time scales Td and T rn for an ex- 
ample case with voltage driving. Assume the parame- 
ters for the oscillator to be L = 200nH and C — 50fF. 
This yields the resonant angular frequency loq = 10GHz 
and characteristic impedance Z b = 2kfL Suppose that 
we are performing the transmission measurement at 
T = Hujo/{4:k b ) = 20mK or T = T%Lo Q /{2k b ) = 40mK, 
corresponding to the red and green curves in Fig. 1101 



We further assume that the amplifier is impedance 
matched such that the amplifier looks to the sample as an 
impedance Zq = Zi,/ 500. This can be accomplished via a 
transformer circuit. Now taking eVb = 0.03?ia;o = 0.2/iV 
and Tjv = 5hu>o/kb — 400mK, we get the measuring 
time T m = 17/cjq — 1.7ns, and the relaxation time 
Td = 93/uiq = 9.3ns. Thus, it should be possible to ex- 
tract the information of the state of the qubit in a single 
measurement. These numbers were calculated by assum- 
ing z = 0.8 and \ (f\p\i) = l/y/2. The ratio T d /T m can be 
further improved for voltage driving if Z b can be made 
larger, Zq smaller, or if the noise temperature Tjv can be 
decreased. Note that the estimate ijTBj) is a pessimistic 
one: the off-resonant mixing time is longer than when the 
driving frequency corresponds exactly to the position of 
the impedance peaks. 
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